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1 Introduction 

The classical compressed sensing theory (see the pioneering work in mm) 
has been gradually generalized and the nonlinear compressed sensing theory 
has attracted more and more attentions inspired by significant real-life applica¬ 
tions such as sparse phase retrieval and sparse pulsation frequency detection 
in Asteroseimology (see |2J and references therein). Among those nonlinear 
measurements, the polynomial structure has been employed in many applica¬ 
tions cases, such as quadratic measurements in sparse signal recovery [25], and 
nonlinear compressed sensing with polynomial measurements [23]. Besides the 
nonlinearity structure carried in many real-life applications, some priori infor¬ 
mation such as the nonnegativity, can be extensively encountered in commu¬ 
nications, DNA microarrays, spectroscopy, tomography, network monitoring, 
and hidden Markov models mmmmm- In this regard, the optimality 
condition such as the KKT condition, which has been widely used in the opti¬ 
mization community, might be a good surrogate for the involved feasible set. 
Under some constraint qualifications, the original problem turns out to find¬ 
ing the sparsest solutions to a special nonlinear complementarity problem with 
polynomial structures. Mathematically, it can be formulated as 

(prp\ min 11*11° 

' s.t. F(x) > 0, x > 0, (x, F(x)) = 0, 

where F = (/i,..., /„) and all /)’s are polynomial functions. 

Apparently, when each _/) in the aforementioned model ( PCP ) is reduced 
to be affine (i.e., a linear function added by some constant), (PCP) is exactly 
to find the sparsest solutions to a linear complementarity problem (LCP for 
short). The sparsest solution to LCP has been studied in [TU1I50] . However, to 
our best knowledge, related work on this topic are very limited, partially but 
essentially due to the complexity caused by the involved non-convex discon¬ 
tinuous objective function, and partially from the nonlinearity generated by 
the underlying complementarity constraint. For the former difficulty, many re¬ 
laxation strategies have been explored by using different surrogates for the £q 
norm, such as the convex t\ norm mmm- the non-convex £ p norm [MM, 
the reweighted i\ norm [BUS] . and so on. A natural but essential question arises: 
is it possible to get an exact solution of the original £q norm minimization 
problem by the relaxation counterpart? If so, what kind of properties should 
the involved data possess? For the linear measurement case, the well-known 
restricted isometry property (RIP for short) was introduced to guarantee the 
desired exactness, which has given a great explanation of the popularity of all 
sorts of random compressed sensing approaches [3M1 - There are some other ex¬ 
act relaxation conditions on the coefficient matrix for linear constraints, such 
as the null space property [88] . the range space property [39], the s-goodness 
property [2T] and so on. Most of these exact relaxation properties are somehow 
not easy to verify. In m, a generalized zT-matrix was introduced. Together 
with the nonnegativity of the right-hand side observation vector, it serves as 
an easy-to-check condition for the desired exact relaxation for the linear sparse 
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optimization. For the nonlinear sparse optimization, such as the sparse LCP, 
the Z-matrix has been employed to guarantee the exact relaxation [50[ . Then 
how about more general polynomial cases? Can we find some Z -type condition 
to ensure an exact solution from the relaxation problem for the original sparse 
polynomial complementarity problem? This is our primary goal in this paper. 

Recently, tensors, as a higher-order generalization of matrices, have been 
extensively studied rnmmm. which is closely related to polynomials in 
terms of those coefficients. This allows us to write a polynomial equation sys¬ 
tem in a condense form with tensors. In this regard, when the constraint func¬ 
tion in the aforementioned problem ( PCP ) takes the form F(x) = H(x) — b 
with some homogeneous polynomial function H(x) and some vector b £ R ra , 
then the feasible set can be reformulated as 

(TCP) x > 0, Ax m ~ 1 - b > 0, (x, Ax m ~ 1 - b) = 0, 

where to — 1 is the degree of H(x), x m ~ 1 is a rank one tensor of order to — 1 
and dimension n with its (zi,..., z m _i)th entry x ^ ■ • ■ Xi m _ 1 , and A is an mth 
order n-dimensional tensor consisting of all the coefficients of H (x) by means 
of H(x ) = Ax m_1 . Here the tensor product Ax m ~ 1 is defined as (Ax rn ~ 1 )i = 

n 

a ii 2 - -i m x i 2 '' ’ x imi f° r all i = 1, ■ ■ ■, n. Similarly, we can define Ax m ~ k 

i • = l 

as 

n 

(*/4.X )ii ...ifc — ^ ^ ...im %ik+i ’ ’ * ^im 5 ^3.5 • • ■ ? — 1, . . . , Tl. 

ifc+i 

_ (i-i) 

The above (TCP) is the so-called tensor complementarity problem which has 
been studied in [313133! ■ In this paper, we will focus on finding the sparsest 
solutions to a tensor complementarity problem which can be modeled as 

(P \ min Iloilo 

1 0) s.t. Ax m ~ l - b > 0, x > 0, (ar. Ax m ~ 1 - b) = 0. 

Mathematically, problem (P 0 ) is generally NP-hard due to the objective func¬ 
tion ||:r||o. Inspired by the scheme of the most popular convex relaxation, we 
could get the a linear surrogate e T x resulting from the nonnegativity con¬ 
straint of x. But the nonlinearity from the tensor complementarity constraints 
cannot be easily handled and the existing exact relaxation conditions are not 
that appropriate since most of them are customized for linear systems. 

In this paper, by employing Z-tensors (see e.g., [131 land the least element 
theory in nonlinear complementarity problems 1341 , we present that if the 
involved b is nonnegative and A is a Z-tensor, then a sparsest solution of 
the tensor complementarity problem can be achieved by solving the following 
polynomial programming problem: 

, „ , min e T x 

^ l ’ s.t. Ax™- 1 


— b = 0, x > 0, 
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where e is the all one vector. In comparison to those existing exact relaxation 
conditions for general nonlinear sparse optimization problems [2], our condi¬ 
tions on the coefficients of the polynomial functions are easy to check. This is 
the main contribution of this paper. 

The rest of the paper is organized as follows. The concepts of Z-tensor and 
Z-function are recalled and some useful properties are presented in Section 
2. The Z-tensor complementarity problem is introduced and discussed in Sec¬ 
tion 3. The exact relaxation theorem is established in Section 4. Concluding 
remarks are drawn in Section 5. 

For convenience of presentation, the following notations will be used through¬ 
out the paper. We use R” and R" to denote the n-dimensional Euclidean space 
and its nonnegative orthant respectively. R" xn is used to denote the space of 
all real nxn matrices. Tm t n is used to stand for the set of all real tensors with 
order m and dimension n. Vectors are denoted by lowercase letters such as x, 
matrices are written as capital letters such as A, and tensors are written as 
calligraphic capital letters such as A. 


2 Z- Tensors and Z-Functions 

As a nonlinear generalization of Z-matrices with non-positive off-diagonal el¬ 
ements, the concept of off-diagonally antitone functions was first introduced 
by Rheinboldt in [23], which further leads to the definition of Z-functions as 
stated in [33]. In [20], Isac has redefined the Z-functions equivalently by means 
of an implication system, which has been also widely used in the community 
of complementarity problems. In this section, the definition of Z-functions and 
some useful properties will be recalled. Particularly, to explore the Z-property 
of homogeneous polynomial functions, the Z-tensor and the partially Z-tensor 
will be introduced and analyzed. 

Definition 1 (Definition 3.2, [20]) A mapping F : R" —> R” is said to 
be a Z-function if for every x, y, z £ R" such that (x, y — z) = 0, we have 
(x,F(y)-F(z))< 0. 

Proposition 1 (Proposition 3.2, [20] ) A Gateaux continuous differentiable 
function F : R” — >• R” is a Z-function if and only if for any x £ R" , VF(x) 
is a Z-matrix. 

Lemma 1 ([20]) If F : R" —> R n is a Z-function, then the following impli¬ 
cation holds: 


x £ R+, y £ R" , (x, y) = 0 => (x, F{y) - F( 0)) < 0. (2.1) 

Moreover, if F{x) = Ax is a linear function, then A is a Z-matrix, which is 
equivalent to the following implication: 


x £ R” , y £ R" , (x, y) = 0 => (x, Ay) < 0. 


(2.2) 
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It is known that for any matrix A , Ax is a Z-function if and only if A is a 
Z matrix, i.e., all of its off-diagonal entries are non-positive. This concept has 
been extended to the higher order tensors as stated below. 

Definition 2 ([37]) Let A = (cii 1 ...i m ) £ 7m,n- zf is called a Z-tensor if all 
its off-diagonal entries are nonpositive, i.e., < 0 when = 0. 

Another concept called the partially Z-tensor is introduced here. 

Definition 3 Let A = (dii...i m ) £ 7 ~ m ,n- We call A a partially Z-tensor if for 
any G [n], < 0 for all i 2 , ■ ■ ■, i m satisfying i x £ {i 2 , ■ ■ ■, i m }- 

Obviously, a Z-tensor is a partially Z-tensor, and both Z- and partially 
Z-tensors of order m = 2 are exactly Z-matrices. Thus, these two concepts 
both can be regarded as extensions of the Z-matrix. Properties on these two 
types of tensors are discussed as follows which will play an essential role in the 
sequel analysis. 

Theorem 1 For any given A G T m ,n, we have 

(i) if A is a partially Z-tensor, then the implication (2.1) holds with F(x) = 
Ax™- 1 ; 

(ii) if A is a Z-tensor, then F{x) = Ax™- 1 is a Z-function. 

Proof, (i) Suppose x, y G R" with (x,y) = 0. Easily we can get 

Xi > 0, t/i > 0,Xiyi = 0,Vi G [n]. (2.3) 


Thus, 

(x,F(y)-F( 0)) 

n 

= ( x,Ay m - 1 ) = J 2 (Ay m ~ 1 ) i 


i= 1 

n n 

= ^ ^ ^ ^ ^722 • • Am y^2 ’ ' ’ Hi- 

7=1 72 ? • • • I'i'rn = 1 

/ 

n 

= £ 

7=1 


^ ^ a ii2 •• -7 m ' ‘ ’ “I - ^ ^ 


&222 -7 m %iVi2 * * ’ V^n 


7 2 ,...,2 m = l 

\7^{7 2 ,...,2 m } 


72 i • • • ?7m — 1 
7G{72,...,7 m } 


n n 


— ^ ^ ' a u 2 ...i m Xiyi 2 • • ■ yi n 

i=[ 

< 0 , 


i=l 


where the last equality follows from (2.3) and the last inequality follows from 
Definition Q5J 
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(ii) Invoking Proposition!]] it suffices to show that for any x £ R", V x ( Ax m_1 ) 
is a Z-matrix. Combining with Lemma Q] we only need to show the implica¬ 
tion (2.2) holds with A := At" 1-1 for any given x £ R”. Let y, z be any two 
nonnegative vectors, and ( y , z) = 0. It yields that 


(■ V , Vi ( Ax m 2 ) z) 


n n 


— ^ ^ Vi ^ ^ 

i=1 i2vi*m = l 

/ 


“(“••• H - *^22 ’ ’ * 1 ^rr 


= E 


< 0 , 


\ 


^ T T ■ • ■ T ) *^2 ' ' ' 1 2/z 




' 2 . ■ ■ ■ j * rri 

\ w 


/ 


where the second equality is from the complementarity of y and z, and the last 
inequality is from the fact that A is a Z-tensor. Thus Fa is a Z-function. □ 


3 Z-Tensor Complementarity Problems 

It is known that a tensor complementarity problem always takes the form 

(TCP{A, b )) X > 0, Ax™- 1 -b> 0, (x, Ax™- 1 - b) = 0, 

which is actually a special nonlinear complementarity problem. When the 
involved tensor A is a Z-tensor, the corresponding ( TCP(A,b )) is called a 
Z-tensor complementarity problem. In this section, we will concentrate on ex¬ 
ploiting the properties of such a special class of tensor complementarity prob¬ 
lems. We start with recalling a nice property possessed by general nonlinear 
complementarity problems with Z-functions. 

Theorem 2 (Ex. 3.7.21, fl7] f Let F : R" — >■ R” be a continuous Z -function. 
Suppose that the following nonlinear complementarity problem 

(. NCP{F)) x > 0,F(x) > 0, (x,F(x)) = 0 

is feasible, i.e., F := {x £ R n : x > 0,F(x) > 0} ^ 0. Then F has a unique 
least element x* which is also a solution to ( NCP(F )). 

Inspired by the relation between Z-tensors and Z-functions, we can easily 
get the following properties for Z-tensor complementarity problems. 

Corollary 1 Let A be a Z-tensor and b £ R ra . Suppose that the tensor com¬ 
plementarity problem ( TCP(A,b )) is feasible, i.e., F := {x £ R n : x > 
0,Ar m_1 — b > 0} 7 ^ 0. Then T has a unique least element x* which is 
also a solution to ( TCP(A,b )). 
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Proof. Theorem [L] tells us that Ax m ~ 1 is a Z-function. Utilizing Proposition 
CD it is easy to verify that Ax m ~ 1 — b is also a ^-function for any b £ R". Thus, 
the desired result follows directly from Theorem [2j □ 

With a nonnegative b and a partially Z-tensor A in (TCP{A, b)), the tensor 
complementarity problem can be equivalent to a multi-linear equation with 
nonnegative constraints. 

Proposition 2 Let A be a partially Z-tensor and b £ R” . The following two 
systems are equivalent: 

(i) x £ R" , Ax™- 1 ~b £ R", (a:, Ax™- 1 - b) = 0; 

(ii) x£Wl, Ax™ -1 — b = 0. 

Proof. Trivially, any solution to system (ii) is a solution to system (i). Let 
y be any solution to system (i). Since A is a partially Z-tensor and b £ R”, 
invoking Theorem [2j it yields that 

0 > (Ay™- 1 ~b, Ay™- 1 ) 

= (Ay™- 1 - b , Ay™- 1 - b) + (Ay™- 1 - b 7 b) 

> \\Ay™~ 1 -b ||i 

This indicates that Ay m_1 — b = 0, which implies that y is a solution to (ii). 

□ 

Note that ■Z’-tensors are partially .Z-tensors. Thus the results in the above 
proposition hold for .Z-tensors. 

Corollary 2 Let A be a Z-tensor and b £ R" . The following two systems are 
equivalent: 

(i) x £ R”, Ax™- 1 — b £ R", (x, Ax™- 1 -b) = 0; 

(ii) x £ R", Ax™- 1 -b = 0. 

Utilizing the aforementioned equivalence, we can characterize the feasibility 
of TCP(A , b) in terms of the consistency of the corresponding nonnegative 
constrained multi-linear equation. Before stating the feasibility, we recall the 
definition of M-tensors, which form an important subclass of Z-tensors. 

Definition 4 ( |3T] ) Let A £ T m ,n be a Z-tensor with A = si — B, where 
I is the identity tensor whose diagonal entries are 1 and others 0, B is a 
nonnegative tensor and s £ R". If s > p(B), then A is called an M- tensor. If 
s > p(B ), then A is called a strong M-tensor. Here p(B) stands for the spectral 
radius of B. 

Proposition 3 If A is a strong M-tensor and b £ R", then ( TCP(A,b )) is 
feasible. 

Proof. Let A = si — B be a strong M-tensor with B > 0 and s > p(B). 
Invoking the equivalence as established in Corollary [2j it suffices to show that 
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that there exists some nonnegative x such that Ax m 1 = b. Let T Sy B,b '■ R+ —> 
R” be the mapping defined as follows: 

T s ,B,b{x) := (s _ 1 £>x m_1 + s^b ) [ ^ ] , Vx e R n , 


where is the vector with its ith component x ^ -1 , for alH = 1 , ..., n. 

Easily, we can find that the required nonnegative solution x is exactly a fixed 
point of this mapping T s ^,b- Besides, since A is a strong M-tensor, applying 
Theorem 3 in ED there always exists a positive z such that Az m 1 > 0. 
Denote 


a := mm 


JLXXJLJLJL , 1 \ 1 

ie[n] (A?™- 1 );’ 


and /3 := max 




Obviously, 


0 KaAz™- 1 <b< fiAz™- 1 . 


/ l \ m -l / l \m-l 

Set v := A I a™- 1 zj and w := A I/3™- 1 z) .Therefore, 
u^=t z = T s # iV < T SyBtb ^zj , 

f3^=r z = T SjB ^ > T Sj b !6 

Note that T s ^,b is an increasing continuous mapping on R" . By employing the 
fixed point theorem in [T] (also see Theorem 3.1 in ED, there exists at least 
one fixed point x of T s ^.b such that 0 < a™- 1 z < x < /3 m ^ 1 z. This completes 
the proof. □ 


4 Exact Relaxation 

Now we are in a position to establish the exact relaxation theorem for the Iq 
norm minimization problem (Po). 

Theorem 3 Let A be a Z-tensor and b £ R” . If the problem (Po) is feasible, 
then (Po) has a solution x* which is also the unique solution to the problem 

(Pi)- 


Proof. Invoking Theorem^ we know that Ax m_1 is a Z-function. Together 
with Proposition [I] it is easy to verify that Ax m_1 — b is also a Z-function. 
Thus, Theorem [2] tells us that there exists a unique least element x* in T as 
defined in Corollary [TJ which is also a solution of the tensor complementarity 
problem. The nonnegativity constraint directly yields that x* is one of the 
sparsest solutions of (Po). Utilizing the equivalence as shown in Corollary [2] 
x* is definitely the unique solution of (Pl) by the fact that it should be the 
least element in {x £ R" : x > 0, Ax m_1 = b}. This completes the proof. □ 
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Corollary 3 Let A be a strong M-tensor and b £ R” . Then problem (Pi) is 
uniquely solvable and the unique solution is also an optimal solution to problem 

(Po). 


Proof. This follows directly from Proposition [3] and Theorem [3] □ 

Some extended result on exact relaxation theory is discussed as follows. For 
any matrix P = ( pij ) £ R raxn , we define a linear operator P.4, : T m ,n —>• T m , n 
as follows: 


( Piii -2 .. .im ^ ^ A — (diii 2 .. -i m ) G "7m,n- 

*=1 

This is also treated as a special tensor-matrix product in the tensor commu¬ 
nity, (see e.g., [23]) Evidently, if P is an invertible matrix, then PT m ,n = T m ,n- 
However, this operator cannot preserve the Z-property for Z-tensors. Note 
that the right-hand side of the multi-linear system Ax m ~ x = b is actually a 
condense form of n homogeneous polynomials of degree m with any row tensor 
Ai formed by the corresponding coefficients of the i th polynomial. If we change 
the order of these equations in the multi-linear system, the solution set will 
not be affected. This observation allows us to generalize the result in Theorem 
[31 For convenience, we use n to denote the set of all tensors of order m and 
dimension n which can be transformed to Z-tensors with some permutation 
matrices, i.e., 

P ^ n := {A G T m ,n '■ PA is a Z-tensor, P is a permutation matrix} 

Corollary 4 Suppose A G P^ n and b £ R” . If the problem (Po) is feasible, 
then (Po) has a unique solution x* which is also the unique solution to the 
problem (P). 


5 Conclusions 

To pursuit the sparsest solutions to a tensor complementarity problem can 
be formulated as an norm minimization with tensor complementarity con¬ 
straints, which is always NP-hard. Based on the properties of iT-tensors, we 
show that one of the sparsest solutions of the Z-tensor complementarity prob¬ 
lem can be achieved in polynomial time by solving a polynomial programming 
problem with a linear objective function, such as the Gauss-Seidel iteration 
method proposed in m- The involved condition on the input of data is indeed 
an exact relaxation condition for the original norm minimization. It is worth 
mentioning that in comparison to other existing exact relaxation conditions 
in the community of sparse optimization or compressed sensing, our proposed 
condition is easy to verify. Including the sparse tensor complementarity prob¬ 
lem as a special case, the topic of sparse optimization with general nonlinear 
complementarity constraints deserves further study. 



10 


Ziyan Luo et al. 


References 

1. Amann, H.: Fixed point equations and nonlinear eigenvalue problems in ordered banach 
spaces. SIAM review 18(4), 620-709 (1976) 

2. Blumensath, T.: Compressed sensing with nonlinear observations and related nonlinear 
optimization problems. Information Theory, IEEE Transactions on 59(6), 3466-3474 
(2013) 

3. Candes, E.J.: The restricted isometry property and its implications for compressed 
sensing. Comptes Rendus Mathematique 346(9), 589-592 (2008) 

4. Candes, E.J., Romberg, J.K., Tao, T.: Stable signal recovery from incomplete and inac¬ 
curate measurements. Communications on pure and applied mathematics 59(8), 1207- 
1223 (2006) 

5. Candes, E.J., Tao, T.: Decoding by linear programming. Information Theory, IEEE 
Transactions on 51(12), 4203-4215 (2005) 

6. Candes, E.J., Wakin, M.B., Boyd, S.P.: Enhancing sparsity by reweighted t\ minimiza¬ 
tion. Journal of Fourier analysis and applications 14(5-6), 877-905 (2008) 

7. Che, M., Qi, L., Wei, Y.: Positive definite tensors to nonlinear complementarity prob¬ 
lems. arXiv preprint arXiv: 1501.02546 (2015) 

8. Chen, X., Ge, D., Wang, Z., Ye, Y.: Complexity of unconstrained I 2 — Ip minimization. 
Mathematical Programming 143(1-2), 371-383 (2014) 

9. Chen, X., Ng, M.K., Zhang, C.: Non-lipschitz-regularization and box constrained model 
for image restoration. Image Processing, IEEE Transactions on 21(12), 4709-4721 
( 2012 ) 

10. Chen, X., Xiang, S.: Sparse solutions of linear complementarity problems. Preprint 
(2014) 

11. Comon, P., Golub, G., Lim, L.H., Mourrain, B.: Symmetric tensors and symmetric 
tensor rank. SIAM Journal on Matrix Analysis and Applications 30(3), 1254-1279 
(2008) 

12. De Lathauwer, L., De Moor, B., Vandewalle, J.: On the best rank-1 and rank- 
(Ri, R 2 , • • • , Rn) approximation of higher-order tensors. SIAM J. Matrix Anal. Appl. 
21(4), 1324-1342 (electronic) (2000) 

13. Ding, W., Qi, L., Wei, Y.: M-tensors and nonsingular M- tensors. Linear Algebra Appl. 
439(10), 3264-3278 (2013) 

14. Ding, W., Wei, Y.: Solving multi-linear systems with A4-tensors. Preprint (2015) 

15. Donoho, D.L.: Compressed sensing. Information Theory, IEEE Transactions on 52(4), 
1289-1306 (2006) 

16. Donoho, D.L., Tanner, J.: Precise undersampling theorems. Proceedings of the IEEE 
98(6), 913-924 (2010) 

17. Facchinei, F., Pang, J.S.: Finite-dimensional variational inequalities and complemen¬ 
tarity problems. Vol. I. Springer Series in Operations Research. Springer-Verlag, New 
York (2003) 

18. Foucart, S., Lai, M.J.: Sparsest solutions of underdetermined linear systems via l q - 
minimization for 0 < q < 1. Appl. Comput. Harmon. Anal. 26(3), 395-407 (2009) 

19. Harmany, Z.T., Marcia, R.F., Willett, R.M.: Spiral out of convexity: Sparsity-regularized 
algorithms for photon-limited imaging. In: IS&T/SPIE Electronic Imaging, pp. 
75,330R—75,330R. International Society for Optics and Photonics (2010) 

20. Isac, G.: Complementarity problems. Springer (1992) 

21. Juditsky, A., Nemirovski, A.: On verifiable sufficient conditions for sparse signal recovery 
via i\ minimization. Mathematical programming 127(1), 57-88 (2011) 

22. Khajehnejad, M.A., Dimakis, A.G., Xu, W., Hassibi, B.: Sparse recovery of nonnegative 
signals with minimal expansion. IEEE Trans. Signal Process. 59(1), 196-208 (2011) 

23. Kolda, T.G., Bader, B.W.: Tensor decompositions and applications. SIAM Rev. 51(3), 
455-500 (2009) 

24. Lauer, F., Ohlsson, H.: Finding sparse solutions of systems of polynomial equations via 
group-sparsity optimization. Journal of Global Optimization pp. 1—31 (2014) 

25. Li, X., Voroninski, V.: Sparse signal recovery from quadratic measurements via convex 
programming. SIAM J. Math. Anal. 45(5), 3019-3033 (2013) 




The Sparsest Solutions to Z-Tensor Complementarity Problems 


11 


26. Liu, Y.F., Dai, Y.H., Luo, Z.Q.: Joint power and admission control via linear program¬ 
ming deflation. IEEE Trans. Signal Process. 61(6), 1327-1338 (2013) 

27. Luo, Z., Qin, L., Kong, L., Xiu, N.: The nonnegative zero-norm minimization under 
generalized Z-matrix measurement. J. Optim. Theory Appl. 160(3), 854-864 (2014) 

28. Qi, L.: Eigenvalues of a real supersymmetric tensor. J. Symbolic Comput. 40(6), 1302- 
1324 (2005) 

29. Rheinboldt, W.C.: On M-functions and their application to nonlinear Gauss-Seidel it¬ 
erations and to network flows. J. Math. Anal. Appl. 32, 274-307 (1970) 

30. Shang, M., Zhang, C., Xiu, N.: Minimal zero norm solutions of linear complementarity 
problems. J. Optim. Theory Appl. 163(3), 795-814 (2014) 

31. Sheikh, M.A., Sarvotham, S., Milenkovic, O., Baraniuk, R.G.: Dna array decoding from 
nonlinear measurements by belief propagation. In: Statistical Signal Processing, 2007. 
SSP’07. IEEE/SP 14th Workshop on, pp. 215-219. IEEE (2007) 

32. Song, Y., Qi, L.: Properties of some classes of structured tensors. Journal of Optimiza¬ 
tion Theory and Applications pp. 1-20 (2014) 

33. Song, Y., Qi, L.: Properties of tensor complementarity problem and some classes of 
structured tensors. arXiv preprint arXiv: 1412.0113 (2014) 

34. Tamir, A.: Minimality and complementarity properties associated with Z-functions and 
M-functions. Math. Programming 7, 17-31 (1974) 

35. Xu, Z., Chang, X., Xu, F., Zhang, H.: regularization: A thresholding representation 
theory and a fast solver. Neural Networks and Learning Systems, IEEE Transactions 
on 23(7), 1013-1027 (2012) 

36. Yang, Y., Yang, Q.: Further results for Perron-Frobenius theorem for nonnegative ten¬ 
sors. SIAM J. Matrix Anal. Appl. 31(5), 2517-2530 (2010) 

37. Zhang, L., Qi, L., Zhou, G.: M-tensors and some applications. SIAM J. Matrix Anal. 
Appl. 35(2), 437-452 (2014) 

38. Zhang, Y.: Theory of compressive sensing via ^i-minimization: a Non-RIP analysis and 
extensions. Journal of the Operations Research Society of China 1(1), 79-105 (2013) 

39. Zhao, Y.B.: RSP-based analysis for sparsest and least ^i-norm solutions to underdeter¬ 
mined linear systems. IEEE Trans. Signal Process. 61(22), 5777-5788 (2013) 




